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the article describes how a normal Shack-Hartmann wave front 
sensor is modified for use in detecting and tracking the focal spot. 
The paper details the analysis and development of the algorithms 
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For direct thrust solar thermal propulsion (STP) or 
terrestrial solar heating using concentrators, proper 
positioning of the solar focal spot is paramount. 

Along with positioning information, a sensor system 
must tolerate the intense heat generated by a solar 
concentrator. The sensor system developed in this 
article satisfies both of the previous requirements. 

The system developed in this article utilizes the 
thruster/absorber as the focal spot sensor with a 
charge coupled device (CCD) camera to image the 
absorber/sensor. Two image processing methods are 
applied to the CCD images to generate the focal 
spot location. Once the focal spot is located on 
an image, positioning information and commands 
are generated to reposition the concentrators. Both 
methods produce the focal spot information needed to 
generate positioning information for the concentrator. 

This article is organized as follows: The STP 
background is presented in Section II. The solar 
sensor problem statement is discussed in Section III. 
Wave front sensors and cylindrical mirrors as applied 
to the problem statement are discussed in Section IV. 
Section V presents the first algorithm developed and 
applied to the model of the concentrator system; 
while Section VI presents the algorithms developed 
specifically for utilizing the wave front sensor. A 
summary and conclusion is given in Section VII. 

II. SOLAR THERMAL PROPULSION BACKGROUND 

The solar concentrator field may be divided into 
two primary areas of interest: terrestrial-based and 
space-based. Both types of concentrator can be either 
imaging or nonimaging and they can be rigid or 
inflatable. Other applications include terrestrial power, 
space power, materials processing, hazardous material 
destruction, building illumination, and space plant 
illumination. 

Solar thermal power as propulsion has its roots in 
the work of Krafft Ericke’s solar powered spaceship 
done in 1956. Ericke’s spaceship was the first design 
that utilized a double solar concentrator connected 
together by a hydrogen powered thruster. 

The major difference between a chemical-type 
thruster or spacecraft and an STP thruster or 
spacecraft is that the STP spacecraft does not carry 
the added weight of the oxidizer as the chemical 
spacecraft. Since the energy for propulsion is 
available to the spacecraft in orbit, in the form of 
sunlight, an STP system would not have to carry 
both propellant and oxidizer to produce energy for 
the upper stage unit. Only an atomic or molecular 
propellant such as hydrogen would be needed on 
the upper stage for thrust. Thus, the STP can have a 
weight advantage over the chemical spacecraft in that 
a lighter molecular weight propellant is used. 
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The STP system heats up hydrogen in the 
thruster/absorber and then expands the hydrogen 
through a nozzle to produce the thrust required by 
the vehicle. Thus the thrust-to-mass (propellant) 
efficiency of the STP system is much better (more 
than two times better) than the equivalent chemical 
upper stage as less nonpayload mass is needed in 
the STP system for the same amount of thrust to be 
developed. [1] Specific impulse or is the rocket 
propulsion term that measures the thrust per unit of 
propellant expended. Thus STP systems have a much 
higher than chemical systems of the same size. 

The definition of is the thrust divided by the 
propellant expelled in the exhaust 


^sp 


gcM 


1 

gcfh 


Jp 

dt 


1 . 

= — -mv. 

gem 


Sc 


( 1 ) 


In (1), Ve is the exhaust velocity. From (1), the 
expression for the power developed by the exhaust 
stream of the thruster is shown in ( 2 ) 


P, = ¥■^^ = \ikglIl. (2) 

Equation (2) can be related to concentrator power 
input to the thruster and the thruster nozzle loss 
mechanism 77 ^ as follows in (3) 


Pe = ( 3 ) 

The solar intensity at Earth is shown in (4) 

IT IT 

/o = 1350^ ±50^. (4) 

The power transferred to the thruster and gas is 
the intensity of the light collected by the concentrators 
minus radiation lost. The power into the thruster is 
then (5) 

= ( 5 ) 

In (5), cr is the Stefan-Boltzmann constant, is 
the intensity from concentrator to absorber A, and T is 
the temperature of the absorber and propellant. In (5), 
is input intensity with effects of the optical path 
through the concentrator acting on it as shown in ( 6 ) 


4 = T-eVo- (6) 

In ( 6 ), Te is the transmission efficiency of the 
concentrator system, A^ is the area of the concentrator 
aperture, and /q is the solar intensity at Earth orbit. 
Defining concentration ratio as = A^/A with A as 
the area of the thruster/absorber face; then from (3), 
(5), and ( 6 ) the power at the exhaust stream is shown 
in (7) 

P, = r,,A{T,C^-aT^). (7) 

Equation (7) is a one-dimensional expression for 
the power into the thruster describing the first-order 
trade-off between concentration ratios and operating 
temperature. The assumptions for (7) are that the 



GOO 


500 


§ <00 
I 

S^300 

c 

c 200 


too 


r (cm) 

Fig. 1. Peak power density versus radius off focus. 
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Fig. 2. Intensity versus solar misalignment. 



thruster has a uniform temperature over A, that the 
backside of A is insulated, and the temperature over 
A is the same for the propellant. Otherwise, (7) would 
become an integral over the flow path. Equation (2) 
and (7) may be used to find thrust F as a function of 
Isp, T, and I^. 

Relating the exhaust power to the misalignment 
angle of the concentrator is a problem solved 
numerically using the concentrator code written by 
Dr. Holmes called offaxis. Eig. 1 shows the result for 
several angles of misalignment. Using values from this 
Eig. 1, Eig. 2 plots the intensity versus misalignment 
angle. Eig. 2 shows that the power density decreases 
considerably at 0.2 deg of misalignment. 

The peak power at 0.2 deg of misalignment is 
approximately 2 dB down and indicates a likely 
cutoff point for this discussion. Therefore, most of 
the work presented in this article pursued a 0.1 deg 
of misalignment target value for maximum power 
transfer to the thruster. 


1308 IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS VOL. 45, NO. 4 OCTOBER 2009 


Authorized licensed use limited to: AFRL. Downloaded on March 25,2010 at 11:17:31 EDT from IEEE Xplore. Restrictions apply. 




















Fig. 3. Solar thermal configuration. 


Misalignment is the difference between the 
position of the Sun and the position of the axis of 
symmetry of the paraboloid-shaped concentrator. 
Misalignment or alignment is really a matter of 
correctly positioning the focused beam from the 
approximately paraboloid concentrator onto the 
absorber so that maximum power is absorbed into 
the propellant. It is not changing the shape of the 
concentrator, but moving struts around in 3-D space to 
position the focal beam onto the thruster. A paraboloid 
reflector is a concave mirror whose cross-sections are 
parabolas or it can be described as a parabola rotated 
about its axis. This type of mirror can concentrate 
light at its focus and provide tremendous heat. 
However, to provide the most heat, the paraboloid 
must be accurately pointed. 

Along with controlling the concentrator during 
thrusting operations, another control requirement is 
that the focal cone should always be positioned on 
the aperture closest to the nozzle for effective heat 
transfer to the propellant. [2] This positioning of the 
focal beam is determined by the misalignment of the 
concentrator with respect to the Sun, and with respect 
to the absorber. 

Since the concentrator for this article is an 
inflatable style, the pointing information derived from 
the sensor designed in this article could be useful for 
providing another input to the inflation control unit 
to increase the accuracy of the paraboloid. This is 
a separate result that was not utilized for the current 
article. 

The current Air Force Research Laboratory 
(AFRL) concept for the concentrators for a solar 
thruster is two off-axis paraboloid concentrators 
connected to a central thruster. See Fig. 3. 

Each concentrator is connected to the central 
thruster with a hexapod or Stewart platform offering 6 
degree of freedom (DOF) control. The 6 DOF needed 
are yaw, pitch, roll, x, y, and z. Coarse alignment 
brings the focal cone into rough alignment with the 


desired aperture on the thruster [2]. Coarse alignment 
is obtained using an on-axis detector, such as a Sun 
tracker or photodiode, which brings the spacecraft 
and concentrators into a general alignment with the 
Sun. Concentrator control should provide optimum 
energy transfer to the solar thruster while protecting 
the concentrator and spacecraft from the dynamic 
behavior of the whole system. 

III. SOLAR SENSOR PROBLEM STATEMENT 

To be useful as a solar sensor especially in space, 
the sensor cannot rely on Earth coordinates for 
tracking information. The sensor also must be able 
to handle the high temperatures involved with a 
solar concentrator that concentrates the sunlight to 
temperatures approaching 3100 K on the absorber. 

The research problem is to: 1) develop a sensor 
that provides information about the location of 
the focus beam from the concentrator that also 
can tolerate the extreme temperatures involved 
in providing solar heating, and 2) develop the 
mathematics of a control system that will focus the 
concentrator in an optimal fashion based on the sensor 
measurements. 

The sensor developed solves the temperature 
issue by using the absorber as a wave front sensor 
for focal spot determination. This sensor overcomes 
the thermal limitations of the other sensors by 
using the reflections of the concentrated light in 
the sensor as the tracking information. Since the 
absorber has to take the high temperatures of the 
concentrated sunlight, it is already capable of handling 
the temperature requirements. The CCD images the 
absorber to determine the location of the focal beam 
based on shape, size, and location of the reflections 
on the absorber. Additionally, the thrust could be 
monitored and used as a feedback mechanism to 
enhance the control of the focal spot. 

IV. WAVEFRONT SENSOR 

The Hartmann wave front sensor was developed by 
Hartmann in 1900 and was used for checking optical 
telescopes for aberrations. It was an array of holes in 
a plate placed in front of the mirror of the telescope 
being checked so that light tubes from the array would 
impinge on the telescope mirror. Two photographic 
plates were used to collect the light information. One 
plate was placed just before the focal spot of the 
telescope and the other one was placed just behind 
the focal spot. The distance between the plates was 
chosen so that the light rays would be separated from 
each other. The path of the light rays was traceable 
by measuring the centroids of the two images. Thus, 
Hartmann was able to determine figures of merit for 
various large telescopes [3]. 
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Fig. 4. Wave front sensor comparison. 


Dr. Roland Shack developed the method of 
replacing the holes of the Hartmann sensor with 
small lenses or lenslets to improve the capability 
of the Hartmann sensor when taking images of 
satellites or stars at the same time as the wave front 
phase error was being determined. The method was 
developed while solving a problem with imaging 
satellites and stars from the Earth. Dr. Aden Mienel’s 
solution was to determine the optical transfer 
function of the atmosphere at the same time as the 
image being taken. The Hartmann sensor was not 
satisfactory as Mienel could not allow the array 
to cover the aperture of the telescope because it 
would eliminate the image of the satellite in favor 
of determining wave front tilt over the imaging of 
the satellites. Mienel tried to use a beam splitter 
to take Hartmann images while taking images 
of the satellites and stars. This method suffered 
from the weak intensity images of the light rays 
in the Hartmann sensor. Shack suggested that 
replacing the apertures in the Hartmann sensor 
with lenses would eliminate the problems that 
Mienel was having with his Hartmann device. The 
lenslets focused spots onto an image plane such 
that the displacement of the spots from the ideal 
position of the spots indicated wave front tilt [3]. 

Thus, the Shack-Hartmann system overcame the 
shortcomings of the Hartmann sensor while increasing 
the sensitivities of the Hartmann device alone. A 
comparison of the two systems is depicted in Fig. 4. 
The Shack-Hartmann sensor is not suitable for use 
in tracking the focal spot of concentrated solar rays 
because of the high temperatures encountered in 
a concentrated sunlight application. As with other 
sensors, lenslets would suffer from the effects of 
concentrated light and from contamination of hot 
absorber vapors. 

The wave front sensor described in this article 
utilizes the tubing of the absorber as mirrors that 
collect information about the light scattered on the 
sensor from the concentrator. In this manner, the tubes 
take the place of the holes or lenslets of the Hartmann 
or Shack-Hartmann wave front sensor. 



Cylindrical mirror (one tube) 




CtKivex i-ciLrmr in licvriidiilal Jir. 

cftTkbf 

ria0.5"Rt 


Fig. 5. Cylindrical mirror. 


The tubes from the conical absorber act like 
cylindrical mirrors (torroidal mirrors) in reflecting 
the sunlight from the concentrator to the CCD. The 
virtual location is important in applying the principles 
from wave front sensing, since the CCD array cannot 
physically be located behind the mirrors or tubing. 

In the current case, the lens system virtually or 
mathematically positions the array behind the absorber 
in a manner similar to the Shack-Hartmann system. 
Each tube would have its own corresponding subset 
of virtual pixels mapped directly onto the CCD array. 
These subsets would not overlap. Fig. 8 shows the 
schematic of the thruster and concentrator with a view 
of the light in the concentrator as seen by the tubes of 
the absorber. 

A. Cylindrical Mirrors 

Since the sensor for this application uses coiled 
tubing as its primary configuration, a brief review 
of cylindrical mirrors is necessary to understand the 
reflections generated by the conical tubing absorber. 
Fig. 5 shows the normal situation encountered with 
the use of a cylindrical mirror. In general a cylindrical 
mirror will tend to compress an object’s reflection 
along the length of the mirror. An object’s reflection 
also tends to get larger as the object moves towards 
the mirror. 

Fig. 5 also shows the reflections from the 
cylindrical mirror along two of its axes. In one axis, 
in the straight on view, the cylindrical mirror acts like 
a flat mirror with the light rays adhering to the rules 
for a flat mirror. That is, incident light is reflected 
such that the reflected ray returns from the mirror as 
if it originated at a point behind the mirror, at equal 
distances to the object; this is also a virtual image. 

Or the reflection occurs such that the reflected ray’s 
angle is equal to the angle of the incident ray. On 
the other axis of the mirror, reflections follow the 
rules of a convex mirror. That is, the mirror surface 
has a virtual focus point and all of the reflected light 
rays appear to come from that focus point behind 
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the mirror surface. Reflections from other portions 
of the cylindrical mirror are more complex than 
the reflections just discussed, however, the overall 
effect is that the reflection in the cylindrical mirror 
is compressed along its main axis and appears larger 
as the viewer moves closer to the mirror. 



B. Flux Tube Calculation 

The flux tube construction shown in Fig. 6 was 
used to determine the ability of the CCD pixels 
to discern changes in light on the concentrator 
due to misalignment via the sensor mirrors. A 
concentrator that is on focus would have light 
filling all of the concentrator as seen by an observer 
located at the focus. The sensor or observer in the 
present case is the absorber/sensor developed in this 
article. 

Since we assume that we don’t lose light from 
the flux tubes, then the pixels should see the same 
intensity from the tube path as it sees looking 
directly at the Sun. Also, we assume that our 
optical system is perfect without any loss of light 
from the flux tubes due to transmission losses or 
reflection losses. The construction presented in 
Fig. 6 illustrates that with perfect optics, the light 
intensity on a pixel via the concentrator mirror path 
is identical to the light intensity that the pixel would 
receive by direct viewing of the Sun. All virtual 
pixels whose solid angle fell within the 30 deg 
half-angle of the cone would have full intensity 
light on them. Therefore, any pixel intensity value 
less than direct Sun intensity and whose solid angle 
fell within the cone of light of the concentrator 
would indicate misalignment of the concentrator. 

If the pixel was at or beyond the angle of the 


X 



Fig. 7. View of figures in cylindrical mirror from camera. 



randam jn-fdcuied posltidn af 
concentrolDr. 

Fig. 8. Light in concentrator. 

cone of the concentrator, that pixel would then be 
considered a border pixel and would not have full 
intensity. 

Fig. 7 shows a schematic of the camera view of 
the image of the symbols and stick figure formed by 
the cylindrical mirrors of the tubing with the CCD 
of the plane of the symbols. While Fig. 8 indicates 
how a CCD matrix would view the reflection of the 
concentrator in the tubing. 
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Fig. 9. Schematic of proposed solution. 


V MODEL SENSOR CALCULATIONS 
A. Methods 

Each image taken could be processed by utilizing 
profile information along diagonal lines representing 
the four quadrants of the circle (along the 45 deg 
angles, say). The areas of maximum intensity would 
be determined along each profile line. The maximums 
should occur roughly where the tubes appear in 
the image as that is where pixel virtual images are 
located. The difference between these areas should 
give an indication of the direction to the focal spot 
(almost a centroiding operation on the maximum areas 
in the image). 

By knowing where the center of the absorber 
is located with respect to the camera (a nontrivial 
assumption as the camera would probably be 
mounted to one of the concentrator’s movable struts), 
the controller should be able to generate x, y, z, 
roll, pitch, and yaw commands for the hexapod 
controller to move the concentrator to a new position 
to provide better focus and thus better heating. 

Fig. 9 shows the schematic of the proposed sensor 
solution. 

For the development of the model presented in this 
article, a copper conical absorber with water-cooling 
if needed is used. Fig. 10 presents the water-cooled 
conical-shaped absorber as shown in the schematic 
in Fig. 9. The fluid to be heated flows through the 
absorber tubing and is heated by the concentrated 
light. 

1312 



Fig. 10. Conical absorber sensor. 


B. Calculations 

Ray tracing was used to determine whether the 
image of the 7 in diameter concentrator model would 
reflect off of the tubes and into the aperture of the 
CCD camera. Fig. 11 illustrates the situation and the 
equation shown could be used in an iteration to find 
the angular extremes that would in fact be imaged in 
the CCD camera. The figure indicates a theta angle 
that would in fact be double the angle with respect to 
the normal vector at the surface of the cylinder. This 
is happening because of the way things were drawn 
in Fig. 11 and the Law of Reflection: the angle of 
incidence is equal to the angle of reflection. A quick 
calculation using just the rays from the extremes of 
the concentrator model using the Law of Reflection 
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Fig. 12. Test apparatus setup. 


showed that the CCD matrix or camera was able to 
view all of the reflection of the concentrator model 
from the tubes. 

The field of view of the CCD camera was 
determined to be 1.4 (495 rows) x 2 (657 columns) 
in at a distance of 1 m from the tubes. 

VI. WAVE FRONT ALGORITHM DEVELOPMENT 

Since the paper [4] presents the hardware being 
used for focal spot tracking in detail, a brief summary 
of the components is now presented for use in this 
article. The schematic of the hardware is shown 
in Fig. 12 and consists of a CCD camera, 1 x 2 m 
elliptical concentrator and a conical-shaped tubular 
absorber. A simulator consisting of a 3 in diameter 
LED taillight was used to represent the Sun’s disk 


at the Earth’s orbit for the 1 x 2 m concentrator 
shown in Fig. 12. The actual value of the solar 
disk at Earth’s orbit is 5 in so the 3 in disk in 
the simulation actually provided a more stringent 
requirement than the operational STP system. The 
Sun simulator was mounted on a tripod that allowed 
various movements of the Sun from on focus to 3-4 
diameters off of on focus for imaging and analysis. 
Putting the Sun simulator at the second focus of the 
concentrator allowed for simulation of the overall 
STP system in a laboratory environment using various 
misalignments off of the second focus point situations 
to be simulated. 

The developed algorithms are based on the 
Shack-Hartmann wave front sensor. The holes and 
lenslets are replaced by the cylindrical mirrors of the 
absorber in the current system. Sunlight reflecting 
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Fig. 13. Flow chart of sequence. Coordinate system after deployment. 


on the absorber tubing is imaged by the camera for 
analysis instead of imaging the array of holes or 
lenslets of the Shack-Hartmann device. The thruster 
tubes serve as high temperature mirrors. The absorber 
takes the heat from the sunlight and conducts it 
to the hydrogen gas for use in providing thrust. 
Differences in location of light and how much light 
is reflected from the tubes are used to determine the 
focal spot location error for use in a feedback control 
system. 

Two or more different tracking algorithms are used 
to locate the focal spot location from the images of 
the tubes. Different algorithms were deemed necessary 
because of the 0.1 deg arbitrary control requirement 
and not having one universal algorithm that covered 
the three regimes of focus coarse focus (sunlight 
more than 5 diameters off focus), intermediate focus 
(approximately 5-3 diameters off focus), and fine 
focus (2-0 diameters off focus). 

A. Method One 

The first method to be discussed is the more 
general method using correlation, while the second 
method is less general, based on area moments 
and area centroids similar to one version of the 
Shack-Hartmann analysis. The flow chart in Fig. 13 
shows the sequence to be implemented for solar 
concentrator focus control. Fig. 14 illustrates the 
situation for the spacecraft just after deployment. The 
coarse focus control system operates on the vector 
positioning the systems so as to minimize this vector. 
The situation then is shown in Fig. 15. The focal spot 
is located in the general area of the sensor but not on 
the optimized spot and that the system is not thrusting. 
Thus, after deployment and coarse alignment, the 
camera is pointed at the absorber and the concentrator 
is positioned facing the Sun. This situation is then 
defined as the initial position for the fine focus control 
algorithms. With the addition of a stored image of 
the absorber in the on-Sun and on focus designed 
configuration, the control algorithms begin. 



Fig. 15. Coordinate system, fine focus control. 

At the start of the control program, the control 
computer has this initial configuration and initial 
image in memory. Also stored in the control computer 
are masks to be used in computing where the focal 
spot is currently located each iteration. The correlation 
or the lack of correlation at the masked areas indicates 
the current location of the spot and provides the 
direction of travel required to reduce the mismatch. 
The overall process for this control algorithm follows 
image correlation or scene matching from computer 
pattern recognition [5-7]. 

From Fig. 15, r^n-or’ P’ *^error 

primary vectors in the development of the algorithms. 

is the vector that determines the location of 
the image in both coordinate systems for the focus 
spot and is the vector to be minimized in the fine 
focus algorithms, is the vector normal to the 
plane at the current image location at p is 
the location of the camera or “eye-position” with 
respect to both coordinate systems. Finally, X(^) 
is defined as the vector representing areas of the 
images that are located in regions of interest. The 
regions of interest in the case of these algorithms 
would be the thruster itself or the individual rings 
of the absorber. The required accuracy of focal spot 
location determines whether the thruster as a whole 
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would be the region of interest or the individual 
rings. 

For the imaging portion of the process, the vector 
X(^) in the preceding paragraph would have to be 
projected onto the CCD camera image plane. The 
method used for this operation is based on pinhole 
optics and projection. The CCD camera utilizes a 
pinhole-like aperture larger than a pixel for its near 
infinite depth of field. A lens is positioned just behind 
the pinhole-like aperture to focus the aperture image 
onto the CCD matrix in the camera. With these optical 
parameters the projection of X(^) onto the CCD 
matrix with coordinates {a,b) is determined using the 
following equations 

a{e) = a * ((X(^) + - p) • e,)/{X{e) + - p) • ^3 

(8) 

b{e) = (3 * ((X(^) + r_ - p) • 2,)/(X(^) + r_ 

(9) 

The parameters a and [3 in (8) and (9) are scaling 
factors determined by the CCD pixel size. The 
e vectors are unit vectors located at p defining a 
coordinate system at p for the camera. 

Now that the image and system coordinates and 
vectors have been defined, the correlation function 
is derived. The image taken by the CCD is denoted 
I(a,b) and is a function of (a,f3) as expected. W(a,f3) 
will represent the image or subimage that we are 
trying to match. For example, W(a,P) would represent 
the ideal on focus image or it could represent only a 
portion of the ideal image such as specific rings of the 
sensor. The L 2 norm is then utilized to calculate the 
minimum distance between W(a,f3) and the current 
image 

/ / / 

E(m,n) = T.Y.VV + m-lj -\-n-l)- 

\i=l ;-l 

( 10 ) 

Where E{m,n) in (10) is the measure of distance at 
coordinates Equation (10) also assumes that W 

is smaller than the image /, but that does not preclude 
using a W that is equal to I in size. Square and 
expand both sides of (10) while using the assumption 
that the energies in W and I are constant and very 
small; the correlation of W and E results as shown 
in (11) 

I j 

E^{m,n) = —2 * + m — 1,7 + n — 1) * 

1=1 j=i 

( 11 ) 

Multiplying (11) by —1/2 generates the correlation 
R(m,n) at (m,n) as shown in (12) 

I j 

R(m,n) = ^^/(/ + m — 1,/ + n — 1) * W(iJ). 

i=l j = l 

( 12 ) 




Tube Blobs" can appear anywhere 
aEong or on the tube. 

Plate "Slobs" appear on the plate 
in the center of the tube. 

Fig. 16. Generic light images on tubes of sensor. 


Finding the maximum value for R(m,n) determines 
the best location in the image that matches W. Finding 
the X,Y location of the maximum value for R(m,n) 
and determining the difference from that location to 
the location of the autocorrelation of the on focus 
image determines which direction to move the image. 
Once the difference between the autocorrelation of the 
on focus image location and the point of maximum 
correlation is found, the control system would move 
the concentrators to reduce this difference. Equation 
(12), correlation, was thought to be sufficient for 
determining the location of the focus. 

B. Method Two 

The second method for determining focal spot 
location is by measuring light area moments in 
the image. This method is derived based on a 
modification of the Shack-Hartmann wave front 
sensor using cylindrical mirrors instead of lenslets 
to register light areas [9]. Images of the concentrator 
can be seen in each coil of the absorber. Whether 
a particular image on the tube, denoted as a “tube 
blob” is bright or dim depends on how close the tube 
is to the focal point. An image area moment using 
the “tube blob” areas is calculated to find the focal 
location. The derivation also follows from a discussion 
of area moments and centroids from mechanics [10]. 

Fig. 16 shows a generic schematic of light blobs 
from an image of the sensor. The shape and intensity 
area of each tube blob is related to the amount of light 
in the concentrator being directed to that area of the 
tube. Calculating the area centroid of the images in 
the tubes determines where the focal spot is located 
on this iteration of the algorithm. Knowing where the 
centroid of the focal spot is located for this iteration 
determines in which direction the control algorithm 
needs to move the concentrator, to move the centroid 
to the center of the sensor. 

Fig. 17 shows the situation that might occur once 
the CCD acquires an image. A coordinate system for 
the image is shown in Fig. 17. The image is then 
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Fig. 17. Light in image from CCD. 



Fig. 18. Centroid calculation. 

processed by first thresholding the image and then 
calculating the area moments. 

Once the image has been thresholded, the area 
moments are calculated and an average or overall 
centroid is calculated as the location of the focal spot 
on the sensor. 

As shown in Fig. 18, the image of the central 
flat plate will have a significant “plate blob” of light 
located at the center of the sensor as the control 
described in this paper assumes a coarse focus is 
obtained before the methods of this paper are used 
for fine focus. If the “plate blob” of light in the center 
of the image is too big for the area centroids to be 
calculated, it can be suppressed in the calculation of 
centroids to obtain locations that are very close to 
focus. 



Fig. 19. Correlation up 1 diameter. 



Fig. 20. Correlation 4 diameter offset. 


C. ImageJ and Image Analysis 

The ImageJ image analysis and processing package 
from the National Institutes of Health [9, 10] was used 
to process the images taken in the experiments. The 
program provides a graphical user interface (GUI) and 
extensive image analysis commands. 

ImageJ performs correlation on images using the 
conjugate multiplication of the Fourier transforms 
of two images to be correlated. Masks were utilized 
to reduce the size of images to include only specific 
tubes of the absorber area as the region of interest. 

Fig. 19 shows the result for the cross-correlation 
between the images for up 1 diameter from on 
focus of the solar simulator and the on focus. In this 
instance, the correlation is smeared and basically 
located at the center of the image indicating that the 
1 diameter misalignment has not been resolved from 
the on focus situation in which the maximum value is 
located in the center of the image. 

Fig. 20 shows a much better situation for 
the correlation technique. This figure shows the 


1316 IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS VOL. 45, NO. 4 OCTOBER 2009 


Authorized licensed use limited to: AFRL. Downloaded on March 25,2010 at 11:17:31 EDT from IEEE Xplore. Restrictions apply. 























Fig. 21. On focus image for area moments. 



Fig. 22. Particles from 21, blobs 1 and 3 are concentrator images 
on tubes. 


correlation between a 4 diameter misalignment image 
and the on focus image. 

ImageJ was also utilized for the analysis outlined 
in method two of this paper. The particular command 
used was the “Analyze Particles” of the Analyze 
menu. To prepare the images for analysis, each 
image had to have a threshold applied to them before 
Analyze Particles would work. 

Fig. 21 presents the situation in the process 
after the image has been cropped and a threshold 
has been applied. The next step in the process 
is shown in Fig. 22 when ImageJ has calculated 
the individual areas based on a minimum size or 
threshold of 100 pixels. The 100 pixel limit came 
from the fact that a 100 pixel area reduced the 
number of areas to work from and was calculated 
to be 3-5 percent of the approximately 2000 pixels 
on a tube blob for a fully lit concentrator on focus. 
Thus we ignored “tube blobs” that were smaller 
than 100 pixels signifying less than 5 percent of 
the concentrator lit. Calculations for the image: 
up0diam_ldiam_pos_x_pt5secs_4_aug_2005 gave an 
average X value of 194.473 and average Y value of 
230.12. The method then uses either the centroid 
of the central “plate blob” of light or the centroid 
of the on focus image to generate the difference 
signal used for concentrator control. For the “tube 
blob” analysis, the difference calculation gives an 


\X\ difference of 64.148 and a |T| difference of 
68.147. By subtracting the focal spot centroids from 
the centroid of the “blob” a direction of movement 
towards focus is generated. When there is no central 
“plate blob” of light in the image, the centroids for X 
and Y from the on focus image have to be used for 
delta calculations. These delta calculations would then 
be used to generate control values for the hexapod 
controller. Looking at Fig. 23, it can be seen that this 
method works down to at least 1.5 in off of focus for 
the solar simulator. That number is associated with 
0.15 deg of offset, which is very close to the arbitrary 
requirement that we control to 0.1 deg. 

VII. SUMMARY AND CONCLUSIONS 

A system was needed to provide focus control 
for solar concentrators used in the high temperature 
environment of terrestrial and space applications. 

The system discussed in this paper performs this 
application superbly. A well-known and used 
method of measuring imperfections in telescopes, 
Shack-Hartmann wave front sensor and in the current 
case, an absorber, was modified and algorithms were 
developed to use in the control portion of the system 
developed. The algorithms were developed using 
ImageJ for processing and analyzing the camera 
images. 
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Fig. 23. Centroid results. 
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The two algorithms developed were the correlation 
and the area moment algorithms. Both methods work 
for locating the focal spot on the sensor down to 
about 2-3 Sun simulator diameters off of alignment. 
This amount of misalignment corresponds to 6-9 in 
off of alignment of the simulator. The 6-9 in of 
offset corresponds to 0.6 to 0.9 deg of angle of 
misalignment. Correlation using only 12 fails to work 
below 1-2 diameters of misalignment. Area moments 
did better below 1-2 diameters of misalignment. Good 
focal spot tracking was seen from 5 diameters of 
misalignment down to 1.5 in of misalignment. The 
1.5 in gives about 0.15 deg of misalignment using 
the 1 X 2 m concentrator, which is very close to our 
original arbitrary requirement. 

The system as developed in this paper solves the 
high temperature problem and gives two algorithms to 
determine focal spot location. 

Future work will include automating, using 
ImageJ, the steps outlined in this paper to control 
either a model of the concentrator and thruster or 
an actual concentrator with hexapod and thruster. 

Also, a combination of both methods should be 
studied to overcome the limitations of the correlation 
method. Finally, data should be taken below 1.5 in 
of misalignment to verify that the method can get to 
the 0.1 deg of misalignment and improve the control 
down to 0.02 deg of misalignment where the thruster 
acceleration could be utilized in the optimization of 
concentrator location using a combination of phase 
only correlation (POC) and area moments or a refined 
area moment method. 
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